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Abstract — We present the results of three-dimensional numerical sim- 
ulations of mass transfer in semi-detached binary systems in which the 
mass-losing star is rotating. The cases of aligned and misaligned non- 
synchronous rotation of the donor star are considered; the resulting flow 
patterns are compared to the synchronous case. The main properties of 
the flow, such as the formation of an circumbinary envelope, the absence 
of a "hot spot" on the edge of the accretion disk, and the formation of a 
shock wave along the flow edge, are qualitatively similar to those obtained 
earlier. For the case of misaligned, non-synchronous rotation, the behav- 
ior of the disk and surrounding matter in established flow regime reflects 
changes in the boundary conditions at the surface of the donor star; in 
other words, a "driven disk" model is realized in the calculations. 



INTRODUCTION 

Since the discovery of non-synchronous rotation (when the rotational velocity of 
one or both components differs from the orbital angular velocity of the system) 
by Schlesinger in the spectral binary 5 Lib in 1909 and in the eclipsing binary 
A Tau in 1910 [1], non-synchronous rotation has been detected in numerous binary 
systems of various types. For a number of Algol-type systems (RZ Set, U Cep, 
RY Per), there is evidence for the rapid rotation of one of the components [2, 3]. 
According to [2, 3], the degree to which the rotation is non-synchronous / = fi^/il 
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is rather high for these systems: / = 6.6 9.0 for RZ Set and / = 10.8 14.5 for 
RY Per. The presence of rotation for one component that differs from the orbital 
rotational velocity has also been detected for several long-period RS CVn stars 
(TZ For, A And, AY Get [4] and a Aur (Capella) [5, 6]). The estimated rotation 
non-synchronicity of the hot star in the a Aur system reaches / ~ 10 12 [7-9]. 

For a number of binary systems, various observational properties (light curves, 
radial velocity curves, etc) also appear to show long-period variations on timescales 
substantially longer than the orbital period. Among the best known binaries of 
this type are the systems Her X-1 (HZ Her) and SS433. In order to explain 
the long-period variations in the Her X-1 system, a "driven disk" model was 
suggested in [10-12], based on the idea that the rotation of the donor star is 
non-synchronous and misaligned. A similar model was also considered later for 
SS433 [13-17]. 

The first attempts to study the role of non-synchronicity in mass transfer were 
made in a ballistic approximation [9, 18]. However, the fact that gas-dynamical 
factors were not taken into account in these calculations casts doubt on the 
reliability of the results. Due to the three-dimensional nature of the problem, a 
correct analysis of the flow pattern in systems with non- synchronous rotation is 
possible only in the framework of 3D gas-dynamical models. Model calculations 
for binary systems with non-synchronously rotating components were presented 
in [19, 20]. Unfortunately, these studies used a simplified formulation of the 
problem, and did not take into account the change of the shape of the donor 
star in the case of non-synchronous rotation [21]; therefore, the results obtained 
require additional verification. In addition, the solution given in [19, 20] was 
limited, since the impact of the circumbinary gas on the flow structure near the 
donor star was not taken into consideration. Calculations made later [22-25] for 
the case of synchronous rotation indicated that the influence of the circumbinary 
gas can substantially change the structure of gaseous flows in the vicinity of L\ 
in a steady-state flow regime. Here, we present results for our solutions to this 
problem formulated in a self-consistent way, devoid of the drawbacks of previous 
models. 

1 THE MODEL 

1.1 Binary system parameters 

In our previous studies [22, 23], we considered the flow morphology in a low- 
mass X-ray binary system with synchronously rotating components in a three- 
dimensional formulation. In order to study the impact of non-synchronous ro- 
tation of the donor star on the flow structure, we now consider a semi-detached 
system with the same parameters as those in [22, 23]. We adopted typical param- 
eters for a low-mass X-ray binary, close to these of X1822-371 [26]. The primary 
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component was assumed to fill its Roche lobe and to have a mass Mi = O.28M0 
and surface temperature T ~ 10^ K; the mass M2 of the secondary, a compact 
object with radius O.O5i?0, was assumed to be 1.4M0. The orbital period of the 
system was Porb = 5''. 56, and the distance between the centers of the components 
was A = 1.97i?0. 

We will use the following terminology. We will call "synchronous" the case 
when all periods — i.e., the rotation periods of both of the stars and of the system 
as a whole — coincide (P+i = P^2 = Porb)', "non-synchronous" rotation is the case 
when the rotation period of the donor star does not coincide with the orbital 
rotation period. For non-synchronous rotation of the donor star, two positions of 
its rotation axis relative to the system are possible: 1) aligned, when the rotation 
axis of the donor star is perpendicular to the orbital plane of the binary system; 
and 2) misaligned, when this axis is inclined to this direction. For all calculations, 
we will assume that the rotation axis of the donor star goes through its center of 
mass and is fixed in the laboratory frame. 

We assumed that the rotation velocity of the donor star in the laboratory 
frame is twice the orbital rotation velocity of the system. Three-dimensional 
calculations were performed for both aligned and misaligned rotation of the donor 
star. 

1.2 Shape of the mass- losing star 

The shape of the mass-losing component in a semi-detached binary system can 
easily be determined for the standard case, when it is assumed that: (i) the orbits 
of the components are circular, (ii) the rotation of the star is synchronized with 
the orbital rotation Q,^ — fl, and (iii) the stars are strongly concentrated, so that 
their gravitational fields can be considered to be the fields of point masses. Under 
these conditions, in the adopted coordinate frame (the X axis is directed along 
the line connected the centers of the stars, the Z axis coincides with the rotation 
axis, and the Y axis forms a right-handed coordinate system whose origin is at the 
center of mass of the donor star), the total potential in the Roche approximation 
can be written 

^ / s GMl GMo 1 /, o\ / s 

^^""^ " " ^ " 2^ (^^ " + ^ ) 

and the shape of the donor star coincides with its Roche lobe; i.e., with the 
equipotential surface passing through the inner Lagrange point Li (see, for ex- 
ample, [27, 28]). When the star reaches the boundary of its Roche lobe, mass 
transfer begins through Li, where the pressure gradient is not balanced by other 
forces. In (1), Mi and M2 are the masses of the components, M — Mi + M2 
is the total mass of the system, Q is the orbital angular velocity of the system, 
fl — (0, 0, Q); Xc — AM2/M; di — ^Jx^ -\-y^ -\- z'^ is the distance to the center of 
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mass of the primary (the donor star), d2 — y{x — A)^ + + z'^ is the distance to 
the center of mass of the accretor, and A is the distance between the components. 

The situation changes when the Roche approximation is not satisfied, so that 
the potential has a shape that differs from (1). Let us consider the general case, 
when the mass-losing component rotates non-synchronously and its rotation is 
misaligned with the orbital rotation: }f Q. Following [29, 30], where the 
potential for the case of non-synchronous rotation is considered, and [21, 31, 32], 
where misaligned rotation is taken into account, we will determine the position 
of the Lagrange points, where the pressure gradient is not compensated by other 
forces. Assuming that fiy^ does not change with time (there is no precession) and 
that the rotation of the star is uniform [21, 33], we can write the equation of 
motion for a test particle in a reference frame rotating with angular velocity fl^, 
in the form (without account for the pressure force): 

, / GMl GM2\ ^ n ^ 

r = - grad — j - Vq - fl^ x [fl^ x r] - 20^ x r , 

where Vq denotes the acceleration of the coordinate origin (the donor star's center 
of gravity) in the laboratory frame: 

ro^nx[nx (-re)] = {n^xc,o,o) 

Vc = (Xc, 0,0) . 

We obtain after simple manipulation 



F=- grad (-^^ - ^Ml \ - O x [O x (-Tc)] - ft^ x [ft^ x r] - 2^1^ x r . 
\ di (12 J 

The positions of the Lagrange points can be determined from the condition 
F = 0. When there exists a potential ^ such that F — — grad^^, these points 
are stationary points of the potential. Using the relations 



-12 X [O X (-rc)] = {-n^Xc, 0, 0) = - grad (O 

-O, X [O, X r] = - grad (-^^^'A^) , 

where A is the distance from the point r to the rotation axis f2*, we find that 
the desired potential ^ has the form 

, , GMl GM2 _o 1^9.9 
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Figure 1: Relative change of the position of the inner Lagrange point for systems 
with non- synchronous rotation Xi^rot / xi^as a function of the mass ratio of the 
components q = M2/M1 for various values of the degree of non-synchronicity of 
the rotation f — Q^/Q. 

and the equation of motion can now be represented 

r = - grad * - 2 fi* x r . 

The potential \1/ includes all forces acting in the rotating reference frame, except 
for the Coriolis force and the pressure gradient. For the special case of aligned 
rotation (O^j. = Q^y = 0) = a;^ + and the potential (to within a constant) 
is described by the expression 

*(-■) = - ^ - io^ ((- - x^f + y') - - o^) (-^ + v') . (3) 

Note that, in the case of aligned, synchronous rotation (O^ = O), the potential 
coincides with (1). 

Using relations (1) and (3), we can determine the change of the position of 
the inner Lagrange point L™* for the case of aligned, non-synchronous rotation 
compared to the standard case. For the sake of convenience, we will write the 
standard Roche potential (1) and the potential for aligned, non-synchronous ro- 
tation (3) in dimensionless form (the distances are divided by A, the potential is 
divided by GMi/A, and q = M2/M1): 
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$(r) = — , , , , - ^ = - -{q + 1) ({x - xcf + y^) , 

Va;2 + 7/2 + ^2 /(a; _ 1)2 + ^2 + ^2 2 ^ ^ 



Figure 1 presents the dependencies of Xirot/xi^ on q = M2/M1 and / = Q^/Q. 
We can see that, when the rotation of the star is slower than the orbital rotation 
(/ < 1), the Roche lobe constructed with account for non-synchronous rotation 
is larger than the standard Roche lobe, with its maximum size achieved when 
f — 0. When the rotation of the star is faster than the orbital rotation, the 
"non-synchronous" Roche lobe is smaller than the standard Roche lobe, with 
Xirot/xLi — > for / — > 00. 

For the case of misahgned rotation, the expression for potential (2) becomes 
very comphcated. We will specify the position of the vector ft^, by two angles: 
the angle between fl^ and the Z axis, and the angle (j) between the X axis and 
the projection of onto the XY plane. The dimensionless potential can then 
be written 

1 



Vx' + y'' + z' ^(x-l)2 + y2 + 
+ q{q+l)x-^{q + l)fA\x,y,z) (4) 



A^(x, y, = - cos^ 0sin^ 1?) - sin^ sin^ 1?) + z^ sin^ ■& 

—xy siv? 1!) sin 20 — xz cos (f) sin 2i3 — yz sin sin 2^? . 

In the case considered, the potential for the misaligned rotation depends on four 
parameters: q, /, i}, and 0. The last parameter changes with time, since the 
vector ft^ rotates with angular velocity — Q in a reference frame rotated with the 
binary system. The inner Lagrange point no longer lies on the line connecting 
the centers of the components. Moreover, according to [21], for some parameter 
values, the potential at the outer Lagrange point L2 can turn out to be lower 
than that at the inner Lagrange point. In this case, the Roche lobe will "open" 
from the side of L2 earher, than it does from the side of Li. However, as this 
occurs for large values of / and (see Fig. 5 from [21]), we will not consider this 
possibility in our study. 
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Let us estimate how much the position of Li deviates from the hne connecting 
the component centers for the case of misahgned rotation. Given that, first, 
ly^rotl < \zLrot \ for not very large i), and, second, the maximum of | Zj^rot I is reached 
when = 0, we wiU estimate the deviation of L™* from the hne connecting the 
component centers using the value \zi^rot\ for (j) — 0. Our analysis of the vertical 
shift of the inner Lagrange point indicates that, for the adopted parameters, 
this deviation is negligible, since it is smaller than the cross section of the flow 
determined using standard models [34, 28]. 

This fact made it possible to simplify the model, and to assume that the 
shape of the star does not change as a function of orbital phase. The shape 
of the donor star is taken to be that of the equipotential surface (4) passing 
through L^°*. Since, for (f) = 7r/2, the inner Lagrange point L™* lies along the 
line connecting the component centers, we used expression (4) for the potential 
with the fixed value = 7r/2 to determine the shape of the donor star. 



1.3 Gasdynamical equations 

The flow of gas in a binary system is described by a system of gas-dynamical 
equations: 



^ + ^ + ^^^ + ^--^f -^S^-P (5) 

dfmj , dpuw , dpyw , d{pw'^ + P) _ _ 5$ 
dt dx dy Bz ~ P~&z 

dpE , dpuh , dpvh , dpwh 9$ 9$ 9$ 



Here, p denotes density; m, v, and w are the x, y, and z components of the 
velocity vector v = {u,v,w); P is the pressure; E = e + ^2 ■ I'^P is the total 
specific energy; h = e + P/ p+ y2-|fpis the total specific enthalpy; and $ is the 
Roche potential. The gas-dynamical equations are written in a reference frame 
that rotates with the binary system, so that the Roche potential (1) is used in 
(5). To close the system of equations (5), we used the equation of state for an 
ideal gas, P = {'j — l)pe. To take into account radiative losses, we took the 
adiabatic index 7 to be 1.01 [35, 36]. The non-synchronicity of the donor-star 
rotation was taken into account when specifying the boundary conditions. 
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1.4 Numerical model 



To solve the above system of equations, we used the TVD scheme of Roe with a 
high approximation order [22, 37, 38]. The system of equations was solved from 
arbitrarily chosen initial conditions to the point when the flow regime became 
steady-state. The calculations for the case of aligned rotation included five orbital 
periods; the solution was stationary over the last three periods. In the case of 
misaligned rotation, it is impossible to achieve a true steady-state regime, due 
to the periodic time dependence of the boundary conditions. Therefore, we took 
the solution to be steady-state when the main features of the flow structure 
repeated with a period equal to that of the boundary conditions. We made 
calculations over six orbital periods, though strict periodicity of the solution was 
already seen for times longer than two orbital periods. The calculation area was 
the parallelepiped [— A...2A] x [—A.. .A] x [—A.. .A]. We used a nonuniform grid 
consisting of 91 x 81 x 55 nodes that were more densely spaced in the accretor 
zone. 

We adopted free-outflow conditions for the matter at the accretor and at the 
outer boundary of the calculation zone. We used (4) to determine the shape of the 
donor star, which coincided with the equipotential surface passing through the 
inner Lagrange point L™* for = 7r/2. The boundary conditions at the surface 
of this star were determined by solving for the decay of the discontinuity between 
the gas parameters (po, Vq, Pq) at the surface of the mass-losing star and in the 
calculation cell that was nearest to the given point of the surface [35] . Note that 
the boundary value of the density does not affect the solution, since the system 
of equations scales in p and P. In the calculations, we chose an arbitrary value 
of Pq] to determine the real densities in a specific system with a known mass- loss 
rate, the calculated densities must simply be increased in accordance with a scale 
determined from the ratio of the real and model densities at the surface of the 
mass-losing component. 

The velocity at the surface of the donor star was specified as 

Vo = 17™* xr + nco, 

where fi™* = fi^ — fi is the angular-velocity vector for the rotation of the donor 
star in the rotating reference frame, cq is the sound speed at the stellar surface, 
and n is the normal vector to the stellar surface. Note that our assumption that 
the star's rotation is uniform implies that the velocity vector at the surface of 
the donor star has a component normal to the surface, which is due, not only 
to the sound speed, but also to the rotation of the star. It is evident that, in 
the equilibrium state for the stellar surface, no non-thermal movement of gas can 
occur normal to the surface. 

Under the assumption of uniform rotation, the star should have the shape of 
an spheroid that extends to L™*. However, the shape calculated using the total 
potential does not yield such a solution, since an additional velocity component 



8 



.CO 



60.0 



40.0 



20.0 



0.0 



-20.0 



-40.0 




4.0 



orb 



Figure 2: Components of velocity vector at the inner Lagrange point for various 
cases: (1) synchronous rotation = cq, Vy = 0, Vz = 0; (2) ahgned, non- 
synchronous rotation = Cq, Vy = x^rotQl"^, = 0; (3) misahgned, non- 
synchronous rotation for -&* — 30° (the angle between the rotation vector in the 
rotating reference frame $7^°* and the Z axis): — Cq, Vy — x j^rotfll"^ cos{'d*) , 
Vz = — a;^rotf2™* sin('(9*) sin(— f2t). The markers correspond to eight times t = 

to + ^sPorb, t = to+ y^Porb, t = to + ^sPorb, ■ ■ ■ , t = to + Porb, -for wMch the 

distributions of gas-dynamical parameters are presented further in Figs. 4c, 7c, 
and 9. 



normal to the surface arises due to the non-spherical shape of the star. This 
implies that the adopted law for the rotation of the star is not self-consistent. The 
problem of determining a self-consistent rotation law for a star in a binary system 
remains unsolved, with the exception of two special cases: synchronous rotation 
and the case when the angular-momentum vector of the donor star is zero in the 
laboratory frame. Simultaneously using in the gas-dynamical model the shape of 
the stellar surface obtained and the gas- velocity boundary conditions determined 
assuming uniform rotation, we obtain a partially self-consistent solution for this 
problem. In the model considered, the gas outflowing from the surface of the 
donor star under the adopted boundary conditions can form new outer layers 
of the donor star, thereby adjusting the self-consistent solution. Unfortunately, 
the presence of additional forces (the Coriolis force and pressure gradient, which 
are not included in the potential) means that we can only approximate a self- 
consistent solution; however, this approach seems optimal for this stage of our 
study. 
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2 CALCULATION RESULTS 



To analyze the impact of the rotation of the donor star on the structure of mass 
flows in the system, we compared models with synchronous rotation (previous 
results, presented in [22-25]), with non-synchronous, aligned rotation, and with 
non-synchronous, misaligned rotation of the mass-losing star. We specified the 
boundary value of the gas velocity at the surface of the donor star using (6); for 
both models with non-synchronous rotation, the rotational velocity of the star 
in the laboratory frame was assumed to be twice the angular rotation velocity 
of the system. For the calculation with misaligned rotation, the stellar rotation 
vector was taken to be inclined to the Z axis by 15° in the laboratory frame (the 
inclination angle in the rotating frame was "d* = 30°). 

The different model calculations differ only in the adopted boundary values 
for the velocity at the surface of the donor star. Figure 2 presents the components 
of the velocity vector at the inner Lagrange point for all cases. In the adopted 
reference frame, the case of synchronous rotation corresponds to the existence of 
only the velocity component V^, while Vy = = 0. In the case of aligned, non- 
synchronous rotation, only Vz = 0, and in the case of misaligned, non-synchronous 
rotation all velocity components are non-zero, with Vz being a periodic function 
of time. The differences in the boundary conditions at the surface of the mass- 
losing star lead, in particular, to changes in the parameters for the matter flowing 
through the vicinity of the inner Lagrange point, as shown in Fig. 3. Figures 3a 
and 3b show surfaces of constant density in the vicinity of the inner Lagrange 
point at the level p = 0.02pLi for the cases of synchronous and aligned non- 
synchronous rotation. The cross section of the YZ plane that forms the edge of 
the plot in Fig. 3 is close to Li (at a distance of 0.066A). The projection of 
the inner Lagrange point onto the YZ cross section is marked by a dot. For the 
case of misaligned, non-synchronous rotation, the adopted boundary conditions, 
and consequently the solution obtained, are time dependent. Figures 3c and 3d 
present surfaces of constant density near the inner Lagrange point at the level 
p = 0.02pirot for the two times t — to+ ^/^Porh and t = to + "^1 iPorbi illustrating 
the change of the flow parameters, and especially the shift of the flow relative 
to the orbital plane. Comparison of Figs. 3a-d shows that the flow parameters 
differ substantially for the different cases, which, in turn, should lead to changes 
of the flow patterns obtained for the different types of donor-star rotation. 

The solutions obtained for synchronous [22-25] and aligned, non-synchronous 
rotation are steady-state. The basic properties of the obtained flow patterns are 
shown in Figs. 4a and 4b, which show surfaces of constant density at the level p = 
0.004pij for both model calculations. For the case of misaligned, non-synchronous 
rotation, we considered the solution to be steady-state when the main features of 
the flow structure repeated with the period of the boundary conditions. Figure 
4c presents surfaces of constant density at the level p = 0.004^^,™* for eight times 
(Fig. 2) covering the total variation period for the boundary conditions. Analysis 
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Figure 3: Surface of constant density in the vicinity of the inner Lagrange point 
at the level p = 0.02^^,^ for the cases of: (a) synchronous rotation; (b) non- 
synchronous, aligned rotation; (c) and (d) non- synchronous, misaligned rotation 
for times t = to + ^/gPorb t = to + ^/ APorh, respectively. The cross section in 
the YZ plane was made at a distance of 0.066A from Li. The thick dot marks 
the projection of the inner Lagrange point onto the YZ cross section. 
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a) 



b) 




Figure 4: Surface of constant density at the level p — 0.004^^^ for the case of (1 ) 
synchronous rotation and (h) non-synchronous, aligned rotation. The thick dot 
marks the position of the accretor. 

of the results presented in Fig. 4c shows that the behavior of the matter in and 
around the disk reflects changes in the injected matter (as well as changes in the 
boundary conditions at the surface of the donor star) . This implies that a driven- 
disk model has been realized in our calculations. An additional illustration of this 
fact is provided by Fig. 5, which shows the temporal dependencies of variations in 
the Z-coordinates of the center of mass of the accretion disk and of the velocity at 
the inner Lagrange point. Analysis of the curves presented in Fig. 5 confirms the 
"driven" character of the solution, and the strict dependence of the flow pattern 
near the accretor on the behavior of the flow in the vicinity of Li . Note that the 
inclination angle of the normal to the disk surface to the Z axis is ±10° in the 
XZ plane and ±20° in the YZ plane. 

Let us consider the details of the flow pattern obtained. For the cases of 
synchronous and aligned, non-synchronous rotation, the main features of the 
steady-state flow pattern obtained are clearly visible in Figs. 4a,b and 6a,b. Our 
results indicate that the rarefied gas of the circumbinary envelope (stream lines 
'a' and '6') has an appreciable impact on the structure of gas flows in the system. 
The gas of the circumbinary envelope interacts with the matter flowing from the 
vicinity of Li and deflects it, leading to a shock- free (tangential) interaction of 
the flow with the outer edge of the forming accretion disk (stream hue 'c'), and, 
consequently, to the absence of a "hot spot" on the disk. At the same time, 
the interaction of the circumbinary-envelope gas (stream line '6') with the flow 
leads to the formation of an extended shock of variable intensity along the edge 
of the flow. The impact of aligned, non-synchronous rotation of the donor star 
is manifested only in a number of quantitative changes, while the overall flow 
pattern remains qualitatively similar. 

As noted above, in the case of non-synchronous rotation, the existence of a 
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Figure 4c: Surface of constant density at the level p = 0.004p^rot for the case of 
non- synchronous misaligned rotation for the eight times indicated in Fig. 2. The 
size of the area displayed for (c) are the same as for (a) and (b). The thick dot 
marks the position of the accretor. -i o 
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Figure 5: Z coordinate of the center of mass of the accretion disk for the case 
of non-synchronous, misaligned rotation as a function of time. The dashed hne 
marks the Z component of the velocity at the inner Lagrange point. 



Vy component for the velocity brings about changes in the parameters of the 
inflowing matter (Fig. 3). Because of the shift of the flow near Li, the flow 
approaches the accretor as it advances, and the disk that forms is substantially 
smaller than in the case of synchronous rotation (Figs. 6a,b). In our calculation 
for non-synchronous rotation, the thicknesses of the flow, the accretion disk, and 
the structures surrounding the disk exceed the disk thickness in the synchronous 
case (Figs. 7a,b; 8a,b). Note that, in addition to forming a shock, the interaction 
of the rarefied gas of the circumbinary envelope (stream line '6' in Fig. 6) with 
the flow also leads to the formation of gaseous "clouds" that precede the front 
edge of the flow outside the disk- formation area (Figs. 4, 8). The appearance of 
these formations is associated with the interaction (collision) of rarefied envelope 
gas outside the equatorial plane with the fiow. The larger thickness of the fiow 
in the case of aligned, non-synchronous rotation brings about the formation of 



Figure 6: (Next page) Lines of constant density and velocity vectors in the equa- 
torial (XY) plane for the case of (a) synchronous rotation; (b) non- synchronous, 
aligned rotation; and (c) non-synchronous, misaligned rotation for time t = 
to + ^s-P- The thick dot marks the position of the accretor. The dashed curves 
depict the Roche equipotentials. The vector in the top right corner corresponds 
to a velocity of 800 km/s. Stream lines for the Rowing matter are denoted 'a', 
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Figure 7: Lines of constant density and velocity vectors in the plane perpendicular 
to the line connecting the centers of the stars (the YZ plane) passing through 
the accretor for the case of (a) synchronous rotation, and (b) non- synchronous, 
aligned rotation. The last isoline corresponds to a density of p — O.OOSpL^. The 
thick dot marks the position of the accretor. The vector in the top right corner 
corresponds to a velocity of 1500 km/s. 



more pronounced "clouds". 

In the calculations for misaligned, non-synchronous rotation, all the basic fea- 
tures of the flow inherent to our previous calculations are maintained. Character- 
istic properties, such as the formation of an circumbinary envelope, the absence of 
a "hot spot" at the edge of the accretion disk, and the formation of a shock wave 
along the edge of the flow, are clearly visible in Figs. 4c and 6c. As noted above, 
a "driven disk" is realized in this case, and there are periodic changes in the 
solution in the steady-state regime (Figs. 4c, 5, 7c). Unlike the case of aligned, 
non-synchronous rotation, the matter in the disk and the surrounding structures 
fluctuate about the equatorial plane, leading to the formation of complicated flow 
structures ("tracks"). The presence of these "tracks" of matter means that the 
disk that forms possesses intermediate linear size and thickness, between those 
for the models with synchronous and aligned, non-synchronous rotation (Figs. 6, 
7, 8). Due to the "driven" character of the solution, the formation of "clouds" 
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Figure 7c: Lines of constant density and velocity vectors in the plane perpen- 
dicular to the line connecting the centers of the stars (the YZ plane) passing 
through the accretor for the case of non-synchronous, misaligned rotation for 
the eight times indicated in Fig. 2. The last isolinc corresponds to a density of 
p = 0.005pirot. The thick dot marks the position of the accretor. The vector in 
the top right corner corresponds to a velocity of 1500 km/s. 
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Figure 8: Surface of constant density near the accretor at the level p — 0.004pL^ 
for the case of (a) synchronous rotation; (b) non-synchronous, aligned rotation; 
and (c) non-synchronous, misaligned rotation for time t = + 'Ys-P- The thick 
dot marks the position of the accretor. The cross section of the surface made by 
the half-planes y = {x > A) and x = A {y > 0) also shown. 
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Figure 9: Surfaces of constant density at the level p = 0.0025p£,rot for the 
case of non-synchronous, misaligned rotation at four times in the vicinity of 
the inner Lagrange point (left) and the accretion disk (right). The sizes 
of the depicted areas are [O..5i?0] x [—ARq.ARq] x [—3Rq..3Rq\ (left) and 
[4Rq..11Rq] X [-3Rq.ARq] X [-2Rq..2Rq] (right). 
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(gaseous formations above and below the equatorial plane) is also periodic (Fig. 
9). Since it is not excluded that these formations can contribute to the total 
emission of the system, together with the emission of the disk and surrounding 
material, the periodicity of their formation should be considered when interpret- 
ing observations. 

CONCLUSIONS 

We have presented the results of three-dimensional numerical modeling of the 
mass transfer in semi-detached binary systems with rotation of the mass-losing 
star. Wc considered the cases of both aligned and misaligned, non-synchronous 
rotation of the donor star, and compared the flow patterns obtained with the 
synchronous case. 

Our analysis indicates that the main features of the flow are qualitatively sim- 
ilar for all the calculations. Study of the flow structure indicates that the rarefled 
gas of the circumbinary envelope has an appreciable impact on the structure of 
gas flows in the system. The circumbinary-envelope gas interacts with the matter 
flowing from the vicinity of Li and deflects it, leading to a shock- free (tangential) 
interaction of the flow with the outer edge of the forming accretion disk, and, 
consequently, to the absence of a "hot spot" on the disk. At the same time, the 
interaction of this gas with the flow forms an extended shock of variable intensity 
along the edge of the flow. 

Unlike the cases of synchronous and aligned, non-synchronous rotation of the 
donor star, in the case of misaligned, non-synchronous rotation, it is impossible 
to achieve a genuine steady-state regime due to the periodic time dependence 
of the boundary conditions. Therefore, in the misaligned, no n- synchronous case, 
the solution was considered to be steady-state when the basic features of the 
flow structure repeated with the period of the boundary conditions. Our analysis 
indicates that the behavior of the disk and surrounding material reflects changes 
in the injected matter (or changes in the boundary conditions at the surface of 
the donor star). This implies that a "driven disk" model has been realized in 
our calculations. The inclination angle of the normal vector to the disk surface 
to the Z axis is ±10° in XZ plane and ±20° in the YZ plane. The calculated 
periodic changes of the shape of the accretion disk and of the surrounding gaseous 
envelope should be observable via the associated changes in the emission of this 
region. 
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